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A simple model of charge ordering is considered. It is shown explicitly that at any deviation 
from half- filling (n 7^ 1/2) the system is unstable with respect to phase separation into charge 
ordered regions with n — 1/2 and metallic regions with smaller electron or hole density. Possible 
structure of this phase-separated state (metallic droplets in a charge-ordered matrix) is discussed. 
The model is extended to account for the strong Hund-rule onsite coupling and the weaker intersite 
i—{ • antiferromagnetic exchange. An analysis of this extended model allows us to determine the magnetic 

C/3 ' structure of the phase-separated state and to reveal the characteristic features of manganites and 

other substances with charge ordering. 

PACS numbers: 71.45.Lr, 75.10.-b, 75.30.Mb, 75.30.Kz 



O 

o 



X 



o 

O i I. INTRODUCTION 

■ The problem of charge ordering (CO) in magnetic oxides attracts an attention of theorists since the discovery of 
^ I Verwey transition in magnetite in the end of thirties Q. An early theoretical description of this phenomenon was 

. given e.g. in ||^. Recently this problem was reexamined in a number of papers in connection with the colossal 
' magnetoresistance in manganites, see e.g. The mechanisms stabilizing the CO state may be different: the 

T-H Coulomb repulsion of charge carriers (the energy minimization requires keeping the carriers as far away as possible, 
similar to the Wigner crystallization), or the electron-lattice interaction leading to the effective repulsion of electrons 
at the nearest neighbor sites. In all cases, charge ordering can arise in systems with the mixed valence if the electron 

■ bandwidth is sufficiently small — large electron kinetic energy stabilizes the homogeneous metallic state. In real 
; materials, in contrast to the Wigner crystallization, the underlying lattice periodicity determines the preferential 

■ types of CO. Thus, in the simplest bipartite lattice, to which belong the colossal magnetoresistance manganites of the 
I . type Ri_xAa;Mn03 (R = La, Pr; A = Ca, Sr), or layered manganites R2-KA:rMn04, R2-22;Ai+2xMn207, the optimum 

' conditions for the formation of the CO state exist for doping x — 1/2. At such value of x the concentrations of Mn'^^ 
^ , and Mn*"*" are equal, and the simple checkerboard arrangement is possible. The most remarkable experimental fact 
here is that even at a; ^ 1/2 (in the underdoped manganites, x < 1/2) only the simplest version of charge ordering is 
experimentally observed with alternating checkerboard structure of occupied and empty sites in the basal plane . 
In other words, this structure corresponds to the doubling of the unit cell, whereas the more complicated structures 
with longer period (or even incommensurate structures) do not actually appear in this case. 
$^ ' Then, the natural question arises: how could we redistribute the extra or missing electrons in the case of arbitrary 
doping level, keeping the superstructure the same as for a: = 1/2 ? To answer this question, the experimentalists intro- 
duced the concept of incipient charge ordered state corresponding to the distortion of long-range CO by microscopic 
metallic clusters 0. In fact, the existence of such a state implies a kind of phase separation. Note that the phase 
separation scenario in manganites is very popular now ||l^. Nowadays, there is a growing evidence suggesting 
that an interplay between the charge ordering and the tendency toward phase separation plays an essential role in 
the physics of materials with colossal magnetoresistance. 

In this paper, we consider a simple model, which allows us to clarify the situation at arbitrary doping. We include in 
this model both the Coulomb repulsion of electrons on the neighboring sites and the magnetic interactions responsible 
for the magnetic ordering of manganites. After demonstrating the instability of the system toward phase separation 
in certain ranges of doping, we also consider the simplest form of the phase separation — the formation of metallic 
droplets in an insulating matrix, estimate parameters of such droplets, and construct the phase diagram illustrating 
the interplay between charge ordering, magnetic ordering, and phase separation. 
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One has to note that the mechanism of the charge ordering considered below (the Coulomb repulsion) is not the 
only one. The electron-lattice interaction can also play an important role, see e.g. p6| . In application to manganites, 
one has to take into account also orbital and magnetic interactions These may be important, in particular, 

to explain the fact that the charge ordering in half-doped perovskite manganites is a checkerboard one only in the 
basal plane, but it is "in-phase" in the c- direction. However, the nature of such charge ordering is not clear yet and 
presents a separate problem: it is not evident that the magnetic interactions responsible for this stacking of afe-planes 
in is indeed the dominant mechanism. Let us also emphasize that the charge ordering in manganites is often 
observed at higher temperatures than the magnetic ordering, and one has to look for a model not relying heavily on 
magnetic interactions. Note that in contrast to magnetic interactions, the Coulomb interaction is one of the important 
factors always present in the systems under consideration. Moreover, it has a universal nature and does not depend 
critically on specific features of a particular system. Consequently, our treatment can be applied also to other systems 
with the charge ordering such as magnetite Fe304 0, cobaltites iQ, nickelates etc. 



II. THE SIMPLEST MODEL FOR CHARGE ORDERING 



Let us consider a simple lattice model for charge ordering: 

H = -t c+cj +V n,nj - /i^n„ (1) 

<i,j> <i,j> i 

where t is the hopping integral, V is the nearest neighbor Coulomb interaction (similar nn repulsion can be also 
obtained via the interaction with the breathing- type optical phonons); /i is the chemical potential, and and Cj 
are one-electron creation and annihilation operators, n.i — cfci. Symbol < i,j > denotes the summation over the 
nearest-neighbor sites. Here, we omit for simplicity spin indices. We also assume the absence of a double occupancy 
in this model due to the strong onsite repulsion between electrons. 

The models of the type (1) with the nn repulsion being responsible for the charge ordering are the most popular ones 
to describe this phenomenon, see e.g. P,pp|j20|] and references therein. Hamiltonian (1) captures the main physical 
effects; if necessary, one can add to it some extra terms, which we will also do in Section V below. 

In the main part of our paper, we will always speak about electrons. However in application to real manganites we 
will mostly have in mind less than half-doped (underdoped) systems of the type Ri_a:Aa;Mn03 with x < 1/2. Thus, 
for a real system one has to substitute holes for our electrons. All the theoretical treatment definitely remains the 
same (from the very beginning we could define operators c and c+ in (1) as the operators of holes); we hope that it 
will not lead to any misunderstanding. 

We consider below the simplest case of square (2D) or cubic (3D) lattices, where for a; = 1/2 the simple two- 
sublattice ordering would take place. As mentioned in the Introduction, this is the case in layered manganites, whereas 
the ordering in 3D perovskite manganites is like this only in the basal plane, the ordering being "in-phase" in the c 
direction. To account for this behavior, apparently a more complicated model would be necessary. 

For the case n ^ 1/2, the model (1) was analyzed in many papers; we follow the treatment of the Ref. [||. As 
mentioned above, the Coulomb repulsion (second term in (1)) stabilizes charge ordering in the form of checkerboard 
arrangement of occupied and empty sites, whereas the first term (band energy) opposes this tendency. At arbitrary 
values of electron density n, we shall at first consider a homogeneous CO solution and use the same ansatz as in 
namely 

= n[l + (-l)V]. (2) 

Such an expression implies the doubling of lattice periodicity, with the local densities ni = n(l+r) and n2 = ri(l — r) 
at the neighboring sites. Note that at n — 1/2 for a general form of electron dispersion without nesting, the CO state 
exists only at sufficiently strong repulsion V > 2t Q. The order parameter is r < 1 for finite V/2t, and the ordering 
in general is not complete, i.e. an average electron density n.^ differs from zero or one even at T = 0. 
We use the same coupled Green function approach as in |^ , which yields 

(E + m)Gi - tkG2 - zVn{l - t)Gi = ^ 

(3) 

[E + /i)G2 - tfeGi - zVn{l + t)G2 - 
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where Gi and G2 are the Fourier transforms of the normal lattice Green functions Gu = {{cicf)) for sites i and 
I belonging respectively to one or different sublattices, z is the number of nearest neighbors and tk is the Fourier 
transform of a hopping matrix element. While deriving (3), we performed a mean-field decoupling and replaced the 
averages (c^Cj) by the onsite densities rij = n[l + (— l)*r]. Solution of (3) leads to the following spectrum: 



E + n^Vnz±^J {VuTzf +tl = Vnz±Uk. (4) 

The spectrum defined by (4) resembles the spectrum of superconductor and hence the first term under the square 
root is analogous to the superconducting gap squared. In other words, we can introduce the charge-ordering gap by 
the formula 



A = VriTZ. 



It depends upon density not only explicitly, but also via the density dependence of r. 
Hence, we get 



A2 + i2. 



(5) 



Note that there is one substantial difference between the spectrum of charge ordered state (5) and superconducting 
state, namely here for n 7^ 1/2 the chemical potential does not enter under the square root in (5) in contrast to the 
spectrum of superconductor where 



Then, we can find the following expressions for the Green functions : 



Gi = 



E+ij,—Vnz—LJk+iO 



G2 



tk 1 
2u)k 27r 



1 



E+fi—Vnz+Wk+iO 
1 



E+n—Vnz—cuk+iO E+fj.—Vnz+uik+iO 



(6) 



where 



47r 



Ak — il , Bk=—{1+— 



47r 



(7) 



After the standard Wick transformation E + iO iE in the expression for Gi , we can find the densities in the 
following form 



ni = n(l + ''")= y 
712 = n(l - t) = / 



1 fpiuJk - M + Vnz) + IH fpi-uJk - fJ' + Vnz) 







[0- 


-) 




-) 


i^kj 



fpioJk - + Vnz) 



^k) 



fpi-ujk - n + Vnz) 



2D.BZ 
d^k 



(8) 



where /f(2/) = 1/ (e^/^ -|- l) is the Fermi distribution function, and Osz is the volume of the first Brillouin zone. 
Summing up and subtracting two equations for m and 722, we get the resulting system of equations for n and /x: 



1 = Vzj ^ [fpi-i^k - M + ynz) - fpiiVk - fJ' + Vnz)] 



J [fri-cok - M + Vnz) + fri^Ok - M + Vnz)] 



rf^k 
20 



(9) 



BZ 



For low temperatures T ^ and n < ^ it is reasonable to assume that 11— Vnz is negative. Hence fpitOk—IJ'+Vnz) = 
and /f(— Wfc — /x + Vnz) = 6{—u)k — A* + Vnz) is the step function. 

It is easy to see that for n = i the system of equations (9) yields identical results for all — A < /i — Vnz < A. From 
this point of view, n = 1/2 is a point of indifferent equilibrium. For infinitely small deviations from n = 1/2, that is, 
for densities n = 1/2 — 0, the chemical potential should be defined as = —A -|- Vz/2 = Vz/2 ■ (1 — r). If we consider 
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a strong coupling case V 2t and assume a constant density of states inside the band, then for a simple cubic lattice 
we have r = 1 — -^^jr, and hence 

where W = 2zt is the bandwidth. Note that for the density n = 1/2 a charge-ordering gap A appears for an arbitrary 
interaction strength V . This is due to the existence of nesting in our simple model. In the weak coupling case V -^21 
and with the perfect nesting, we have A ~ VFexp { — ^} and r is exponentially small. For zV ^ W or accordingly 
for V ^ 2t: Aw Vz/2 and r — > 1. As mentioned above, for a general form of electron dispersion without nesting the 
charge ordering exists only if the interaction strength V exceeds certain critical value of the order of bandwidth W 
1^. Further on, we restrict ourselves only to the physically more instructive case of strong coupling V 3> 2t. 

Now let us consider the case n = l/2 — S, where <5 ^ 1 is a deviation from density n — 1/2. For this case /i = /i((5, r), 
and we have two coupled equations for /x and r. As a result. 

Correspondingly, the energy of a charge ordered state is as follows 

Eco{S)=Eco{0)~—S-^S' + O{S'), (12) 

where EcoiO) energy precisely for density n = 1/2 and 1^170(0)1 ^ W. At the same time, the 

charge-ordering gap A is given by: 



2 



2W 



(13) 



Actually, the dependence of the chemical potential /i and the total energy E on 6, Eqs. (11), (12), stems from this 
linear decrease of energy gap A with deviation from half-filling. 

For n > 1/2 the energy of charge ordered state starts to increase rapidly due to the large contribution from Coulomb 
repulsion (the upper Verwey band is partially filled for n > 1/2). Contrary to the case n = 1/2, for n > 1/2 each 
extra electron put into the checkerboard CO state necessarily has occupied nearest neighbor sites, increasing the total 
energy on As a result, we have for \6\ = n — 1/2 > 

Eco{5) = Eco{0) + [Vz - — j 1^1 - -^S' + 0{S'). (14) 
Accordingly, the chemical potential has the form 

M^) = V^.-^-^l^l+o(^^). (15) 
It undergoes a jump equal to Vz for r — s- 1. Note that the gap A is symmetric for n > 1/2 and is given by 



Vz 
"2" 



2W^ 



We could make the whole picture symmetric with respect to n = 1/2 by shifting all the one-electron energy levels and 
the chemical potential by Vz/2, i.e., defining ji' = ji — Vz/2. This change of variables, of course, would not modify 
our conclusions. 



III. PHASE SEPARATION 



The most remarkable implication of (11)-(15) is that the compressibility k of the homogeneous charge ordered 
system is negative for densities different from 1/2, 
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loc^=-^ = — ^-^<0 (16) 
K dn dS dS^ 3Vz ' 

where d = 1/2 — n. This is the manifestation of the tendency toward phase separation characteristic of the charge 
ordered system with S ^ 0. The presence of a kink in the Eco{S) (cf. Eqs. (12), (14)) imphes that one of the states, 
into which the system might separate, would correspond to the checkerboard CO state with n = 1/2, whereas the 
other would have a certain density n' smaller or larger than 1/2. This conclusion resembles that of (see also 
1^,0), although the detailed physical mechanism is different. The possibility of phase separation in the model (1) 
away from half-filling was also reported earlier in ||lj] for the infinite-dimensional case. Below we concentrate our 
attention on the situation with n < 1/2 (underdoped manganites); the case n > 1/2 apparently has certain special 
properties — the existence of stripe phases etc. |l^, the detailed origin of which is not yet clear. 

It easy to understand the physics of phase separation in our case. As follows from (13), the CO gap decreases linearly 
with the deviation from the half-filling. Correspondingly, the energy of the homogeneous CO state rapidly increases, 
and it is more favorable to "extract" extra holes from the CO state, putting them into one part of the sample, while 
creating the "pure" checkerboard CO state in the other part of it. The energy loss due to such redistribution of holes 
is overcompensated by the gain provided by the better charge ordering. 

However, the hole-rich regions would not be completely "empty," like pores (clusters of vacancies) in crystals: we 
can gain an extra energy by "dissolving" in them a certain amount of electrons. By doing this we decrease the band 
energy of the electrons due to their delocalization. Thus, this second phase would be a metallic one. The simplest 
state of this kind is a homogeneous metal with the electron concentration Um- This concentration, as well as the 
relative volume of the metallic and CO phases, can be easily calculated by minimizing the total energy of the system. 
The energy of the metallic part of the sample is given by 

= -~tzn,^ + ct{n^f'^ -f V{n,^f (17) 

where c is some constant. 

Minimizing (17) with respect to n^, we find the equilibrium electron density in the metallic phase. For the strong 
coupling V > zt, get 

n™o = tz/2V (18) 

Thus, in this simple treatment, the system with rimo < n < 1/2 would undergo phase separation into the CO phase 
with n = 1/2 and the metallic phase with n — rimo- Relative volumes Vm and vco of these phases for arbitrary n can 
be found from the Maxwell construction: 

Vm/wCO = (1/2 -")/("- "mo), (19) 

from which we find that the metallic phase occupies the part Vm of the total volume v given by the relationship 

v,n/v^{l/2-n)/{l/2-n^o), (20) 
The metallic phase would occupy the whole sample when the total electron density n is less than Uma- 



IV. AN EXAMPLE: THE PHASE SEPARATED STATE WITH METALLIC DROPLETS 

As we argued above, the system with the short-range repulsion (1) is unstable toward phase separation for n close 
but different from 1/2. The long-range Coulomb forces would however prevent the full phase separation into large 
regions containing all extra holes and the pure n = 1/2 charge ordered region. We can avoid this energy loss by 
forming, instead of one big metallic phase with many electrons, the finite metallic clusters with smaller number of 
them. The limiting case would be a set of spherical droplets, each containing one electron. This state is similar to 
magnetic polarons ( "ferrous" ) considered in the problem of phase separation in doped magnetic insulators mi^,|ll| . 

We present below the estimate for the characteristic parameters of these droplets. The main aim of this treatment is 
to demonstrate that the state constructed in such a way will have the energy lower than the energy of the homogeneous 
state, even if we treat these droplets rather crudely and do not optimize all their properties. In particular, we will 
make the simplest assumption that the droplets have sharp boundaries and that the charge ordered state outside these 
droplets is not modified in their vicinity. This state can be treated as a variational one: if we optimize the structure 
of the droplet boundary, its energy would only decrease. 
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The energy (per unit volume) of the droplet state with the concentration rid of droplets can be written in total 
analogy with the ferron energy in the double exchange model (see [MUiy| ) . This yields 



E, 



droplet 



-trid z 



i?2 



6Vz 



R 

ndT-n- I — 
3 \ a 



(21) 



Here, a is the lattice constant and R is the droplet radius. The first term in (21) corresponds to the gain in kinetic 
energy of electron delocalization inside the metallic droplets, and the second term describes the charge ordering energy 
in the remaining insulating part of the sample. 

Actually, one should include the surface energy contribution to the total energy of the droplet. The surface energy 
should be of the order of W^R^ /V . For large droplets, this contribution is small compared to the term oc R^ in (21); 
it would also be reduced for a "soft" droplet boundary. One can show that even in the worst case of a small droplet 
(of the order of a few lattice constants) with the sharp boundary, this contribution would not exceed about 20 percent 
of the bulk contribution. That is why we will ignore this term below. 

Minimization of the energy in (21) with respect to R gives 



R 



1/5 



(22) 



The critical concentration Udc corresponds to the configuration where metallic droplets start to overlap, i.e. where 
the volume of the CO phase ( the second term in (21)) tends to zero. Hence, 



6 fa 



oc 



V) 



(23) 



By comparing (12) with (21), (22), we see that for the deviations from the half-filling 0<(5<(5c = l/2 — Udc the 
energy of the phase separated state is always lower than the energy of the homogeneous charge ordered state. The 
energy of the droplet state (21) with the radius given by (22) is also lower than the energy of the fully phase separated 
state, obtained by the Maxwell construction from the homogeneous metallic state (17). Correspondingly, the critical 
concentration Udc (23) is larger than Umo (18). There is no contradiction here: the droplet state, which we constructed 
has electrons confined in spheres of radius i?, and even when these droplets start to overlap at n = Udc, occupying the 
whole sample, the electrons in this state, by construction, are still confined within their own spheres and avoid each 
other. In other words, in our droplet state certain degree of charge-ordering correlations is still present, decreasing 
the repulsion and hence the total energy. 

Thus, the energy of a phase separated state with the droplets corresponds to the global minima of the energy for 
all < (5 < 5c- This justifies our conclusion about phase separation into charge ordered state with n = 1/2 and a 
metallic state with small spherical droplets. 

The situation here resembles that of partially filled strongly interacting Hubbard model, with the CO state corre- 
sponding to an antiferromagnetic state of the latter and with the nearest-neighbor interaction V playing the role of the 
Hubbard's U . In both cases, the kinetic energy of doped carriers tends to destroy this "antiferro" or charge ordering, 
first "spoiling" it in their vicinity and finally leading to the formation of the metallic state (Nagaoka ferromagnetism). 
In the Hubbard model, we also face the situation with phase separation at a small enough doping [2l| ]. 

Note also that for n > 1/2 the compressibility of the charge- ordered state is again negative i — = — < 
and has the same value as for the case n < 1/2. As a result, it is more favorable again to create a phase-separated 
state for these densities. However, as it was already mentioned, the nature of the second phase with n > 1/2 is not 
quite clear at present, and therefore we do not consider this case here. 



V. AN EXTENDED MODEL 



Now we can extend the model discussed in the previous sections by taking into account the essential magnetic 
interactions. In manganites, besides the conduction electrons in eg bands, there exist also practically localized t2g 
electrons, which we now include to our consideration. The corresponding Hamiltonian has the form 

H = -t ^ cf^Cjcr + V ^ UiTij ~ JH^^Siffi + J ^ S,,Sj -^^ni, (24) 

<i,j>,a <i,j> i <i,j> » 
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In comparison to (1), the additional terms here correspond to the strong Hund-rule onsite coupHng Ju between 
locahzed spins S and the spins of conduction electrons a, and a relatively weak Heisenberg antiferromagnetic (AFM) 
exchange J between neighboring local spins. In real manganites, the AFM ordering of the CE type in the CO phase 
is determined not only by the exchange of localized t2g electrons but to a large extent by the charge- and orbitally- 
ordered eg electrons themselves. For simplicity, we ignore this factor here and assume the superexchange interaction 
to be the same both in the CO and in the metallic phases. 
It is physically reasonable to consider this model in the limit 

JhS >V >W > JS'^. 

In the absence of the Coulomb term, this is exactly the conventional double exchange model (see e.g. Note that 

the absence of doubly occupied sites in (20) is guaranteed by the large Hund's term. It also favors the metallicity in the 
system, since the effective bandwidth in our problem depends upon the magnetic order. Therefore, the estimate for the 
critical concentration changes here in comparison to (23). Similar to ||T^ the metallic droplets will be ferromagnetic 
(FM) due to the double exchange. The energy of one such droplet has the form 



E 



-t z 



7rV\ 










6Vz 














-zJS^ 















(25) 



The last two terms in (25) describe respectively the loss in the energy of the Heisenberg AFM exchange inside the 
FM metallic droplets and the gain of this energy in the AFM insulating part of the sample. The minimization with 
respect to the droplet radius (as in (21)) yields 



R 



2 \ -1/5 



(26) 



Note that at t/V ^ JS^/t, formula (26) gives just the same estimate for the radius of FM metallic droplet 
(R/a) ~ {t/JS^y/^ as in |Jl|. 

In the opposite limit when {t/V) ^ JS^/t, we reproduce the same result {R/a) ^ {V/t)'^^^ as in (22). Finally, 
critical concentration ric is estimated as follows 



t 

rin oc I — 



3/5 



(27) 



As a result, taking into account also the tendency to the phase separation at very small values of n [p|- pl]Jl^ we 
come to the following phase diagram for the extended model (cf . : 

1. At < n < I ) , it corresponds to the phase separation into a FM metal with n — n' > Q embedded in 

the AFM insulating matrix {n = 0). To minimize the Coulomb energy, it may be again favorable to split this 
metallic region into droplets with the concentration n' and an average radius given by Eq. (26) with t/V = 0, 
each containing one electron and kept apart from one another. 



2. 



At < n < (^-f + < 1/2, the system is a FM metal. Of course, we need a window of parameters 

to satisfy the inequality in the right-hand side. In actual manganites where t/V ~ 1/2 -^1/3 and J/t ~ 1/3, 
these conditions upon n are not necessarily satisfied. Experiments suggest that this window is present for 
Lai_j;Caa;Mn03, but it is definitely absent for Pri_a;Caa;Mn03 |11[; 



3/5 



3. Finally, at 



3/5 



< n < 



we have the phase separation in the form of FM metallic droplets inside the 



AFM charge ordered matrix. 



This phase diagram is in a good qualitative agreement with many available experimental results for real manganites 
p2|- p5[ |, in particular with the observation of the small-scale phase separation close to 0.5 doping p^. Note also that 
our phase diagram has certain similarities with the phase diagram obtained in 
ferromagnetism in doped excitonic insulators. 



29 1 for the problem of spontaneous 
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VI. CONCLUSION 



Summarizing, we have shown that the narrow-band system, which has the checkerboard charge ordering at n 1/2 
(corresponding to the doping x = 0.5) is unstable toward phase separation away from half-fiUing (n ^ 1/2). It 
separates into the regions with the ideal CO (n = 1/2) and the other regions, in which extra electrons or holes are 
trapped. The simplest form of these metallic regions could be spherical metallic droplets embedded into the CO 
insulating matrix. Simple considerations allow us to estimate the size of these droplets and the critical concentration, 
or doping Xc — 1/2 ~ at which the metallic phase would occupy the whole sample and the CO phase would 
disappear. The account of the magnetic interactions does not change these conclusions but somewhat modifies the 
characteristic parameters of the metallic droplets. 

The long-range Coulomb interaction may also modify the results, but we do not expect any qualitative changes. 
For the realistic values of parameters, the size of metallic droplets is still microscopic (about 10-30 A), and the excess 
charge in them will be rather small. 

The obtained picture corresponds rather well to the known properties of 3D and layered manganites close to (less 
than) half doping, x < 1/2. Percolation picture of transport properties emerging from this treatment is confirmed 
by the results reported in |l7|, |l5| , p2|j24| -^ ; moreover the coexistence of ferromagnetic reflections and those of the CE 
type magnetic structure typical of the CO state at x = 0.5 were observed by the neutron scattering ||2^. Thus, the 
general behavior of underdoped manganites {x < 0.5) is in a good qualitative agreement with our results. 

Our treatment leads to the same tendency to phase separation (instability of the homogeneous CO phase) also for 
overdoped regime, x > 0.5. What would be the second phase in this case, is not yet clear. Therefore we did not 
concentrate our attention on such a situation. 

Our treatment is applicable also to other systems with the charge ordering, such as cobaltites ||l8| and nicke- 
lates It would be interesting to study them for charge carrier concentrations different from the commensurate 
"checkerboard" one 

A number of important problems still remain unsolved (the origin of the "in-phase" ordering in perovskite man- 
ganites in c-direction, the detailed description of inhomogeneous states in overdoped regime n > 1/2, the behavior at 
finite temperatures). Nevertheless, in spite of the introduced simplifications, our model seems to capture the essential 
physics underlying the interplay between phase separation and charge ordering in transition metal oxides. 
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